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1. Introduction* 

In  the  context  cf  a general  form  of  the  multivariate  linear  model,  we  consider 
a sequence  { X^  = (X^, . . . ,X  .) 1 , i>^  1}  of  independent  random  vectors  (i.rv)  with 
continuous  cumulative  distribution  functions  (cdfj 


Fi(x)  = P{  £ x}  = F (x-  a-  3ci),i  >l,x  c Rp,  p >1, 


(1.1) 


where  a = (a,,...,  a )'  and  3 = ((  B-J).  , _ _ (q  > 1)  are 

unknown  parameters, 'and  = (c^, . . . ,c  ^)  1 , i>_  1 are  vectors  of  known  regress- 
ion constants. Tiie  form  of  the  cdf  F is  not  assumed  to  be  specified. By  reference 
to  the  usual  canonical  reduction  of  the  multivariate  linear  hypothesis  [viz., 
Anderson (1958), Ch. 8 J , we  consider  the  following.  We  partition 


3 = ( §i,  32),  q1+q2  = q,  0 and  q2  > 0 , (1.2) 

pxq  pxqx  pxq2 
and  consider  the  null  hypothesis 

Hq:  82  = 0 against  the  alternative  H^:  32  i 0 . (1.3) 

Note  that  whenever  q^  is  >_  1,  both  and  11^  are  composite  hypotheses . 

For  the  particular  case  of  q^  = 0 i.e.,  for  8=0,  the  X^  are  identically 
distributed  (i.d.)  under  Hq  and  a class  of  conditionally  as  well  as  asymptotic- 
ally distribution- free  rank  order  tests  has  already  been  studied  by  Puri  and 
Sen(1969).  For  q^  >_  1,  under  Hq*.  82  = 0 , the  X^  are  no  longer  i.i.d.,  and  this 
invalidates  the  approach  of  the  above  mentioned  paper.  In  fact,  in  such  a case, 
genuinely  distribution- free  rank  order  tests  may  not  generally  exist.  In  the 
univariate  case  (i.e.,  for  p=l  ),  this  difficulty  has  been  circumvented  in 
some  specific  problems  by  Sen (1969)  and  Puri  and  Sen (1973)  by  using  suitable 
aligned  rank  order  tests  where  the  alignment  is  based  on  rank  order  estimates 
of  the  nuisance  parameters. This  approach  is  systematically  explored  and  developed 
here  for  the  general  multivariate  linear  model  , and,  basically,  the  theory  of 
rank  order  estimators  of  regression  parameters,  developed  in  Sen  and  Puri (1969) 
and  Jureckova’(1971) , is  employed  here  for  the  estimation  cf  the  nuisance 
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parameters  (i.e. , 3^)  and  incorporated  in  the  construction  of  suitable  aligned 

rank  order  statistics  on  which  tests  for  = 0 are  based. In  the  univariate 

parametric  case,  the  classical  likelihood  ratio  test  possesses  some  (asymptotic) 
optimality  properties.  The  picture  is  somewhat  different  in  the  multivariate 
general  linear  models  where  the  Jikelihcod  ratio  tests  may  not  perform  uniformly 
better  than  others,  even  asymptotically.  However,  under  quite  general  regularity 
conditions,  the  V.'ald-optimality  (viz.,  best  average  power  and  the  most  stringency) 
applies  to  the  likelihood  ratio  type  tests.  In  the  nonparametric  multivariate 
problems  too,  a variety  of  tests  is  available  in  the  literature  [viz.,  Puri 
and  Sen (1971) 3*  and  some  of  these  can  be  adapted  for  the  general  linear  model 
under  consideration. Among  other  possibilities,  we  'consider  here  a class  of 
rank  order  tests  having  some  analogy  with  the  likelihood  ratio  type  tests. For 
such  rank  tests,  it  is  possible  to  develop  distribution  theory  in  some  closed 
forms  comparable  to  those  of  the  likelihood  ratio  tests  and  permitting  an 
asymptotic  comparison  of  these  tests  for  local  alternatives. 

The  proposed  rank  order  tests  for  Hq  are  considered  in  Section  3 following 
the  preliminary  notions  and  basic  assumptions  in  Section  2.  Section  4 deals  with 
asymptotic  comparisons  of  parametric  (mostly,  normal  theory  and  likelihood  ratio) 
and  rank  order  tests.  Asymptotic  optimality  of  the  proposed  aligned  rank  tests 


is  also  considered  in  this  section.  The  last  section  deals  with  a special  case 
of  (1.3),  namely,  testing  the  hypothesis  of  parallelism  of  several  (multiple) 
regression  surfaces,  which  turns  out  to  be  the  multivariate  multiparameter 
analogue  of  Sen (1969). 

2.  Preliminary  notions 

Let  R. . = Zn  . u(X..  - X.  ) (where  u(t)  = 1 or  0 according  as  t is  > or  <0) 

ji  s=l  v ji  jsy  v v b — 

be  the  rank  of  X..  among  the  X.  , s=l , . . . ,n,for  i=l, . . . ,n, j=l, . . . ,p.  Since  F 
J i j a 

is  continuous,  ties  among  the  observations  may  be  neglected,  in  probability. 

For  each  j(-l p)  and  n(^  f),  we  consider  a set  of  real  valued  .scores 

a^  (1) , . . . ,a^  (n) , generated  by  a score  function  fj(u),0  <u  <1,  in  either 
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of  the  following  ways: 

a^(i)  = 4>j(i/(n+l))  or  (uni) » i=1»  • . • ,n;  j=l,. . . ,p,  (2.1) 

where  4 j (u)  is  assumed  to  be  square  integrable  inside  (0,1),  and  <...  < 
U are  the  ordered  random  variables  of  a sample  of  size  n from  the  rect- 
angular  distribution  over  (0,1).  Our  proposed  testing  procedure  is  based 


on  the  following  type  of  rank  order  statistics: 

S = ((S  .,));  S .,=  Ln  .(q.  - q )a^(R..) 

-n  n.jk'-'’  n,jk  1=1'-  ki  lor  n x ji' 

where  c^  = n * c^,  for  k=l,...,q;  j=l,...,p.  We  denote  by 


(2-2) 


— -1  „n 

c = n L . . c ■ 
~n  i=l  ~i 


Following  Ha'jek  (1968)  and  Hooffding  (1973),  wo  assume 
that  for  every  j(=l,...,p)  , 


(u)  = tp^(1)  (u)  - g^2)  (u) 


(2.3) 


where  (u)  , s = 1,2  is  absolutely  continuous  and  ncn- 

decreasing  in  up  (0,1)  and 


a.  / » j 

f |ca  j (u)  ! (u(l-u)  T'-du  < a,  ; s = l,2;  j = l P . (2.4) 

•0  3 


Denote 


: . = j g f u) du  , j = 1, . . . ,p  , 
J > ^ 3 


(2.5) 


00  OP 

xjj'(p)  *J  hi  <ptH  ,X,,V  <F[  j']  <y) »«[  jj.j  (x.y)  -~fy  (2.6) 

where  ^ (x)  and  Fr  (x ,y)  are  the  marginal  cdfs  of  the 
jth  and  the  (j,j')th  components  respectively. 


Assume 
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A(F)  - ((  A^,(l:)))  is  positive  definite (p. d. ) and  finite. 


(2.7) 


We  also  denote  by 

5n  = «Cn,kk'»  -'^-iCSi  * £n)(  Si'  ^ 

((  zi=i Ccki  - Cj^Kcjj.n  - ck,n)  )) 
and  assume  that 

is  p.d.  for  every  n _>  nQ  , 

and  there  exists  a p.d.  and  finite  matrix  C,  such  that 


(2.8) 


(2.9) 


-b 


n C -*•  C as  n + “ . 


(2.10) 


Further,  we  assume  that 


c.  = cV  - c , for  i = l,...,n, 

~i  ~i  ~i  * » » > 


(2.11) 


.(s) 


where  for  each  k (=l,...,q)  and  s(  =1,2),  is  non-decreasing  in  i. 


This  assumption  is  a slightly  simplified  version  of  a parallel  assumption 
made  by  Jureckova(1971) . For  q = 1,  (2.11)  is  not  necessary. 

Finally,  we  assume  that  for  every  e >0,  there  exists  an  integer  n^ 


= n^(e),  such  that  for  n _>  n^. 


max 


n^C  > e { , ' I c,  . 

n,kk  1 <i  <n  1 ki 


c, I },  k 1, . . . ,n. 


"kn1 


(2.12) 


Regarding  the  c.d.f.  F,  we  assume  that  for  each  j (=1, . . . ,p) , the 


marginal  c.d.f.  F^  has  an  absolutely  continuous  probability  density 


function  (p.d.f.)  f^  with  a finite  Fisher  information 


Ij  = Ufj-jj)  = C ((d/dx)log  f .j(x)}“dF  j (x),  j=l,...,p.  (2.13) 

As  has  been  mentioned  in  Section  1,  our  testing  procedure  rests  on 
some  aligned  rank  order  statistics.  To  explain  the  alignment  procedure, 
we  need  the  following  notations. 


Let  B = ((b.,)). 


be  a real  matrix,  and  we  write 


jk  j 1 » • • • > p > k 1 1 • • • t c{ 

B'  = ((  bx bp  ))  where  bj  = (bj  1 b^) , for  j=l Let  then 
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X.(B)  = Xi-Bci  = (Xu(b1),...,Xpi(bp))'  for  i=l,...,n; 


(2.14) 


Rj i = RjiCbj)  = zJ=1u(X:jiCbj)-  XJs(bj))  , i=l . . ,n; j=l, . . . ,p,  (2.15) 

so  that  R. .(B)  is  the  rank  of  X..(b.)  among  X.  (b .) ,s=l, ... ,n,  for  i=l,...,n. 
ji  - ji  ~ j js 

Now,  replace  the  in  (2.2)  by  Rj^(bj)  for  i=l,...,n,  j=l,...,p  and 
denote  the  corresponding  matrix  of  rank  order  statistics  by 


5.®  * «- Sn,jk(S;»h-l PiW q- 


(2.16) 


Note  that  Sn(B)  iii  (2.16),  viewed  as  a function  of  the  pq  elements  in  B, 
generates  a pq- dimensional  stochastic  process  (on  RF^  ).  We  shall  make  use 
of  the  same  in  the  next  section  to  introduce  the  proposed  aligned  rank 
order  statistics. 

3.  The  proposed  aligned  rank  order  tests 


As  in  (1.2),  we  introduce  the  following  partitionments: 

B = ( , §2  ) ’ Si  = ($1(1) » $1(2)^  * ^=l»’*’»n* 

pxq  pxq^  pxq^  qxl  q-^xl  q?xl 
Then,  under  Hq  in  (1.5),  we  have 

Fi(x)  = F(  x - a - ),  i=l,...,n. 


(3.1) 


(3.2) 


First,  we  proceed  to  estimate  the  nuisance  parameter  (matrix)  8^  for 
the  model  (3.2).  For  this,  we  consider  the  pxq1  matrix 


~n(l)^l’^  " ^ Sn,jkQj  ^j=l p;k=l,...,q1  * 


where 


(bj1'))  = (bj1»...,bj  ,0,...,0),  j=l,...,p. 


(3.3) 


(3.4) 


Now,  under  (3.2),  Sn^(  3^,0)  has  the  same  distribution  as  of  Sn^(0,0) 
under  the  hypothesis  that  3=0;  for  the  latter  case,  we  may  use  the  results 
of  Puri  and  Sen(1969)  and  obtain  the  following;  (a)  under  (3.2),  ( 3^,0) 

has  expectation  0 and  dispersion  matrix 


? 


A(P)  M C ,,,,  where  C = 
~ J ~n(ll)  ~n 


~n(ll)  ’ ~n(12) 

~n(23)  ~n  (22) 


(3.5) 


is  of  order  q^xq^  , j.,j=l,2  and  & stands  for  the  Kronecker  product, 
and  (b)  as  n 00  , 


£<  n'h  W > -►  XCxtt,  (?,  }<«»  £(ii)  ). 


(3.6) 


where  is  the  upper  q-^xq^  principal  minor  of  C,  defined  by  (2.10). 

To  estimate  8j,  we  adopt  the  alignment  procedure  studied  in  detail  by 
Sen  and  Puri (1969)  and  Jureckova(1971) , and  define 

?n  " {?1  : £)=1  Zk=l  I * } • (3 

Then,  our  estimator  of  3^  , under  the  model  (3.2),  is  given  by 

A 

6.  = centre  of  gravity  of  D . (3 

— _L  j 1 1 ~I1 


(3.8) 


By  arguments  parallel  to  those  of  Jure£kova'  (1971),  it 


follows  that 


where 


sup  ||8-  - 0 ■ j|  5 0 , as  n 


BCD 
— 1 — n 


~1  -l.n 


,n^i.„  -5P>  -JC„(o.  r(F)  ®c  > 


(3.9) 


(3.10) 


T(F)  = ( ( T , . / (F)  ) ) - (U..,(F)/A.A.,)) 
01  31  1 1 


(3.11) 


Aj  = J (d/dx)cp..  (Ff  (x)  JdFj.  (X)  , j = l p.  (3.12) 


8 


Under  the  model  (5.2),  2,  is  a translation- invariant,  robust,  cons is - 

~ ■* 

tent  and  asymptotically  normally  distributed  estimator  of  8-j . U’e  use  the 
same  for  our  alignment  process  and  consider  the  following  aligned  rank  order 


statistics: 


~n(2)  ^Sn,  jk^j=l, . . . .p;k=q1+l, . . . ,q 


(3.13) 


where 


Sn,jk  Zi=l(cki  ' ckn^an^^ji^’  ql+1",,,q’  (3-14) 


and  the  aligned-ranks  are  defined  by 

R..  = R. . ( 8.  , 0)  for  i=l,.-..,n;  i=l,...,p. 

ji  ~l,n’  ~ J * 

To  introduce  the  proposed  test  statistics,  we  first  define 


(3.15) 


Mn  " ^mjj',n^j,j'=l,...,p 

= ((  (n-1)"1  Z"=1  4j)(Rji)anj,)(Rj'i)  " ^ny)  (3-16) 


where 


= n'1  Z"=1  a^(i)  , for  j=l,...,p. 


(3.17) 


Also,  replacing  R^  by  R^  and  R^,^  by  R^ , ^ for  i=l,,..,n  and  j,j'=l,...,p 
in  (3.16),  we  denote  the  corresponding  matrix  by 


M = ((m..,  ))..,.. 

n ^ JJ'  ,n^j,j'=l,...,p 


Let  then 


* -1 
?n  = ~n(22)  " ~n(21)  ?n(ll)  ?n(12)  * 

A A ^ 

G = M S C (of  order  pq,  x pq,  ) 
—11  —11  —11  **  *- 


(3.18) 


(3.19) 


(3.20) 


~n  ^ Sn,jk  Sn,j'k'-*-*j,j,=l,...,p;k,k,=q1+l (3---0 

where  lln  is  also  of  the  order  pq2  x pq,.  Tlien , our  proposed  test  statistic 


9 


C = Tr(  H G 
- rP  vP 


-1 

n 


(3.22) 


= riy  ylJ  yH  yn  c c 

j=l"j,=l  “k=q1+li'k'=q^+l  n,jkn,j'k'  "n 


/'ii ' *kk' 
nrJ  c 


n 


where 


M'1  = ((m^'j).  . and  C*'1  = ((c*kk')), 

~n  n ,p  ~n  n ■'•'k.k'=q^+l, . . ,q 


(3.23) 


The  analogy  of  to  the  classical  Lawley- Hotel ling  trace  criterion  for 
the  multivariate  analysis  of  variance  (MANOVA)  problem  can  readily  be 
identified.  Whereas  the  latter  is  based  on  the  least  squares  estimators  of 
the  parameters  involved  in  the  model,  ourjCn  is  based  on  the  corresponding 
rank  order  estimators  of  the  nuisance  parameters.  Further,  the  asymptotic 
equivalence  of  the  Lawley-Hotelling  trace  and  the  likelihood-ratio  criteria 
(under  the  normal  theory  model)  is  well-known  [vi;.,  .Anderson  (1959)  ,C’n.  8j  . 

In  the  same  spirit,  we  could  have  proposed  an  alternative  test-statistic 

Xn  “ I l i/l !h„  * s„i l • <3-24) 

where  | | A | | stands  for  the  determinant  of  the  matrix  A.  By  using  Theorem  3.1 
of  Jure"ckova  (1971)  and  proceeding  as  in  the  proof  of  Theorem  3.5  of  Puri, 

Sen  and  Gokhale  (1970),  it  can  be  shown  that  *£n  and  -21og£°  are  asympto- 
tically equivalent  under  the  null  hypothesis  and  for  local  alternatives  too. 

As  such,  in  the  sequel,  we  shall  be  mainly  concerned  with  the  statistic  £ . 

In  the  reminder  of  this  section,  we  show  that  Under  H^:  B-,  = 0,  when 
the  assumptions  of  Section  2 are  met,  £n  has  asymptotically  a chi-square 
distribution  with  pq-,  degrees  of  freedom.  This  provides  an  ADF  (asymptotically 
distribution- free)  test  for  Hq. 

Lemma  3.1.  Under  11^:  g,  = 0 and  the  assumptions  of  Section  2, 


2.  • 1 p 

" 5n  — 

^ A'V)  8 C 1 , as  n 

-*•  00  » 

(3.25) 

* 

- r r c • C = 

If (2 1 ) .(11)  ~ 

/~(H) 

-(12)  \ 

S = £(22) 

.(21) 

J 

£(22)  ) 

| . (5.26) 

where 
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Proof.  By  virtue  of  (2.10),  to  prove  (3.25),  it  suffices  to  show  that 


M 4 /\  ( F ) , as  n 
— n 


(3.27) 


Also  since  m . . , - m . 

33  »n  33 


(j),,,  -(j).2 


= (n-1)  [ y fa  -1  (i)  - a 

]3.n  li41  n n 


-*  x .(F)  = > • by  (2.1)  and  some  routine  commutation^,  we 

33  '33 

need  only  to  show  that  for  every  j / j ' , 


m..,  5 \ ../(F)  when  II  holds  . (3.28) 

33  .n  33  0 

By  assumption  (2.3),  (see  also  Ha'jek  (1968),  section  5) 
for  every  e > 0 , there  exists  a decomposition 

cp_.  (u)  = cp  ! ^ (u)  + cp  (u)  - cp ^ (u)  , 0 < u < 1 , (3.29) 

j = 1 p 

where  cpf^  is  a polynomial,  cpf2^  and  cp  . are  non- 
3 3 3 

decreasing,  and 


3 1 


y [ [cp!k)  (u)  l2du  < e X ■ ■ > 1 < j < p • (3.30) 

J 3 33  '■ 


k“2  0 


Using  (3.29)  we  decompose  m. into  9 terms.  Using  the 

33  / n 

Cauchy -Schwarz  inequality  for  the  eight  terms  for  which  at 
least  one  factor  is  non  polynomial  along  with  (3.30),  it  follows 


that  to  prove  (3.28),  it  suffices  to  take  cp  . - cr  . , 1 < j < p . 

3 3 

Since  the  cpf^  are  absolutely  continuous  and  are  polynomials, 

for  them,  the  corresponding  m..,  can  be  written  as 

3 3 > n 

/“  rm  *i1)0Ls(x))  .,(y))dfi  *.,(x,y)  + o(l)  where  fl  is  the  sample 

■°°  “°°  J nj  j iij  J 


■ 


n 


c.d.f.  for  the  aligned  observations  on  the  jth  variate,  j=l,...,p  and  "njji 
is  the  bivariate  sajnple  c.d.f.  (for  tiie  jth  and  j'th  variates)  for  these 
aligned  observations.  For  the  aligned  vectors  X^(gj,0),  i=l,...,n,  the  sample 
c.d.f. 's  arc  denoted  by  11^  and  H^,  > for  j (/  j')  = l,...,p.  Then,  by  (2.14), 
(3.10)  and  the  continuity  of  the  parent  c.d.f,,  it  follows  that  as  n -*»  , 
sup{|  iinj  j,  (-x,y)  ~ Hnj  i » (x,y)  I X,y  e R2  } •>  Oa.s.,  for  every  j^j  '=1, ...  ,p. 
Also,  note  that.  the.  $ . are  bounded,  continuous  functions.  So,  first  replacing 

A is 

Hn  by  Hn,  11^  by  and  then  using  Theorem  4.1  of  Puri  and  Sen(1969),  the 
desired  result  follows.  In  fact,  it  can  be  shown  that  (5.27)  holds  a.s. 

Lemma  3.2.  Under  = 0 and  the  assumptions  of  Section  2,  as_  n -*•  °°  , 

5.(2)  - 5.(1) ( h-°j  * AJ  h,n  - 5p?n(12))  S'  (3-3« 


where  defining  the  Aj  by  (3.12), 


(3.32) 


A = Diag(  A1,...,  A ) . 

The  proof  follows  as  a direct  multivariate  extension  of  Theorem  3.1 
of  Jure^kova'  (1971),  and  hence,  the  details  are  omitted. 

By  noting  that  cf . Jureckova’”  (1971)  S^  ^ ( £2^,0)  = 0^(1),  the  following 
lemma  also  follows  from  Tlieorem  3.1  of  Jureckova(1971) . 

Lemma  5.  3. Under  the  hypothesis  of  Lemma  3. 2, as  n 00  , 

n ' {§n(l)(  §1’2)  ‘ ~l,n  " 2P  2n(ll)  } 2 ’ (3.33) 

From  (3.31)  and  (3.35),  we  arrive  at  the  following 
Lemma  3.4.  Under  the  hypothesis  of  Lemma  5.2  , as  n °°  , 

n ^ §n(2)  ‘ 5n(2)(  §1»2>  + ^(l)^!’^  ~n(ll)  2n(12)^  ^ 2 ' 
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mean  0 and  dispersion  matrix 

A0:)  « Cn  , (3 

it  follows  that  under  Hq  in  (1.3), 

c>C(n  §n(2)t?l*?)  ' 5n(l)^l*^  Cn(ll)  ~n(12)  ^ 

— > ^pxq2(  A(F)  Q (22 j " £(2l)?(ll)S(l2)}  )> 

Hence,  by  (3.34)  and  (3.36),  we  obtain  that  under  Hq  in  (1.3),  as  n 

1,(2)  ) — ’ Kxqd  8 V-  » 


(3.35) 


(3.36) 


(3.37) 


From  (3.37),  Lemma  3.1  and  the  asymptotic  distribution  theory  of  quadratic 

forms  associated  with  asymptotically  multinomial  vectors,  it  follows  that 

under  the  hypothesis  of  Lemma  3.3,  £n  has  asymptotically  a chi  square 

distribution  with  pq0  degrees  of  freedom. 

Thus,  the  proposed  ADF  test  of  size  ct(0  < a <1)  is  as  follows: 

2 

Let  Xpq  a be  the  upper  100ct%  point  of  the  chi  square  c.d.f.  with  pq-, 
degrees  of  freedom.  Then,  the  null  hypothesis  11^  in  (1.3)  is  accepted  or 
rejected  according  as  £n  is  < or  ^ . 

4.  Asymptotic  comparison  with  parametric  counterparts 

We  confine  ourselves  to  local  alternatives  for  which  the  power  of  the  proposed 
tests  are  away  from  0 and  1.  We  consider  a sequence  {K^}  of  Pitman- type 
alternative  hypotheses,  viz.,  for  some  (fixed)  non-null  y, 

Kn  : h = ?2(n)  = n 2 12'12  =((:Yjk))j=l,..,p;k=q1+l,...,q  (4,1) 

for  which  an  asymptotic  power  function  can  be  traced  and  compared  with  the 

parallel  function  for  some  parametric  tests  for  the  same  problem. 

For  the  normal  theory  model  (where  the  underlying  c.d.f.  F is  assumed  to 

be  multivariate  normal),  classical  parametric  ( likelihood  ratio,  Lawlcy- 

Notelling  trace  or  the  largest  characteristic  root  criterion  of  S.N.Roy  - 

sec  Chapter  8 of  Anderson (1959))  tests  are  all  based  on  the  least  squares 

estimators.  Sen  and  Puri  (1970)  have  studied  the  asymptotic  properties  of 
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the  likelihood  ratio  (as  well  as  the  Law ley -Motel ling  Trace)  statistic  when 
the  underlying  F is  not  necessarily  normal.  It  follows  that  if  !■  possesses 
a finite  and  positive  definite  dispersion  matrix  1(F)  = ((cov(Xj  ■ ,X^. , ^) )) 

= ((  cJjjt  ()')))»  then,  under  in  (1.3),  (=  - 2] og- likelihood  ratio  criterion 

when  F is  assumed  to  be  normal)  has  asymptotically  a chi-square  disdtribution 
with  pq?  degrees  of  freedom.  Also,  under  (K  } in  (4.1),  L has  asymptotically 
a non-central  chi  square  distribution  with  pq7  degrees  of  freedom  and  non- 
centrality parameter 


\ = Trace ( T (1(F)  0 C )_1  ), 

Li  — — 


where 


~ ^ Yjk  Yj'k,^j,j'  = l,...,p;k,k,=q1+l,...,q  ’ 

and  the  y ^ are  defined  in  (4.1). 

t. 

Consider  now  a sequence  of  alternatives  {K^}  , specified  by 

Kn : 3 = ( 0,  n"'5  Y2  ) > Y2  definod  in  (4-1).  (4.4) 

Then,  ( (§p ,?) » §n(2)^l'^  * im<^er  > has  the  same  joint  distribution 

as  that  of  under  . Noting  this  fact,  using  the  results  of  Puri  and  Sen 
(1969)  and  our  lemmas  in  Section  3,  it  follows  by  seme  routine  computations 


that  under!  Kn),  as  n °°  , 


l**  i„(2)  > X™  ( • ««  0 £*  >• 


(4.5) 


From  Lemma  3.1  and  (4.5),  we  conclude  that  under  {K^}  , has  asymptotically 
a non-central  chi  square  distribution  with  pq7  degrees  of  freedom  and  non- 


centrality parameter 


A,  = Trace ( T ( T(F)  8 C*)'1  ), 


(4.6) 


where  T(F)  is  defined  by  (3.11). 


From  (4.2)  and  (4.6),  we  conclude  that  the  Pitman  asymptotic  relative 


efficiency  (ARL)  of  with  respect  to  L^  is  given  by 
CH,L  = hC.'\  = tr(T(VfF)a:*)'b/tr(T(I(F)8C*)'1) 


(4.7) 


which  depends  on  r , F and  C*.  If  F is  multinoraal  and  for<C  we  use  the  normal 


[ 
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scores  (i.e.,’all  the  4>j  being  tl'*e  inverse  of  a standard  normal  c.d.f.),  then, 


it  can  easily  be  checked  that  T(F)  = Z(F),  and  hence,  (4.7)  reduces  to  l,i.e., 
the  aligned  rank  order  normal  scores  test  and  the  likelihood- ratio  test  are 
asymptotically  power- equivalent  for  normal  F and  local  alternatives  in  (4.1). 
However,  in  general,  for  arbitiary  F,  ^ is  bounded  as  follows: 


chp(  E(F)T'1(F)  ) < e^,  < ch2(  Z (1)T  1 (F) ) 


(4.8) 


where  ch. (A)  stands  for  the  jth  (largest)  root  of  A for  j _>  1.  The  bounds  in 


(4.8)  may  be  studied  as  in  Sen  and  Puri  (1967)  or  Puri  and  Sen(1969),  and 
hence,  the  details  are  omitted.  For  testing  simple  hypotheses  in  multivariate 
linear  models,  Puri  and  Sen  (1969)  have  stuied  (in  their  Theorem  6.2)  the 
optimality  of  rank  order  tests  for  local  alternatives.  In  passing,  we  may 
reinark  that  under  the  same  set  of  regularity  conditions  as  in  Theorem  6.2 
of  Puri  and  Sen(1969),  £n  has  asymptotically  the  best  average  power  with 
respect  to  suitable  surfaces  in  the  parameter  space  (of  Y2),  it  has  also 
asymptotically  the  best  constant  power  on  such  surfaces  and,  finally,  it  is 
an  asymptotically  most  stringent  test. 

5.  ADF  tests  for  parallelism  of  several  regression  surfaces  . 


As  a multivariate  generalization  of  the  univariate  problem  treated  in  Sen 

• 00 


(1969),  we  consider  here  the  following.  Let  , k=l,.,.,n^  be  n^  independent 
rv's  with  continuous  c.d.f. ’s 

F{k)(x)  = P{  X.M  £x)  = F(x  -ak  - 3R  cP°),  (5.1) 

for  i=l,...,n,  , k=l,...,sQ>  2).  We  desire  to  test  the  null  hypothesis 

Hq:  8|  = ...  = 8g  = 6 (unknown)  (5.2) 

where  and  B are  treated  as  nuisance  parameters.  If  we  let  3^  = 

8l  + £k  , k-l,...,s  (so  that  p*  = 0),  q = st  ( where  each  of  the  c^  in 
(5.1)  is  a t-vector),  then,  we  arc  in  a position  to  use  the  theory  dcveloved 
in  Sections  3 and  4.  Therefore,  without  giving  the  details  of  derivations,  we 
present  the  main  results  in  this  case  as  follows. 


1 


«;  t 1 


as  in  (2.2J 


r 


Based  on  the  observations  in  the-  kth  sample,  we  define  S^ 
fki 

and  S';  ; (B)  as  in  (2. 14)- (2 . 16) . Let  then  n = n,+...+n  and 

- IS 

S (B)  = if  . (B)  for  B e Rpt. 

-n  ~ k=l  ~n^ 


(5.3) 


Under  Hg,  we  estimate  the  common  B as  follows:  as  in  (3.7)- (3.8),  we  let 
?n  { ? : *j=l  ill  ^n.jm^j^  = ”^imum}  ; (5.4) 


Bn  = centre  of  gravity  of 


Let  then 


i(k)  = sfk)  ( BJ  , k=l,...,s; 


A ' 


H 00  « rr  £00  £00  „ 

~nk  U nk’jr  nk»j ’ r’ j » j ’=1, - • • ,p; x ,r’ =1, . . . , t 

C(k)  = E.k  ( Jk)  - c )(  cPO  - c )’ 

~n,_  i=l v ~i  ~n,/ v ~i  ~n  J 


M„  - M'1^  J) ; (5.9) 

= M 8 CW 


(5.5) 

(5.6) 

(5.7) 

(5.8) 


9n 


k 


~n  ~n, 


for  k=l,... ,s, 


(5.10) 


r(j) 


where  the  c and  aU7  are  the  averages  of  the  regression  vectors  and  the 

scores,  and  are  defined  as  in  Section  2;  rP^  is  the  rank  of  xPO  -E*  ,g  . cP^ 

ji  ji  r=l  n,jr  lr 

among  the  n^  aligned  observations  on  the  jth  variate  in  the  kth  sample,  for 
i=l,...,n^,  j=l,...,p;  k=l,...,s.  Then,  the  aligned  rank  order  test-statistic  is 

4 ■ L1™  9'1  )•  o.ii) 

Under  Hg  in  (5.2),  has  asymptotically  a chi-square  distribution  with  p(s-l)t 


degrees  of  freedom,  while  under  the  sequence  of  alternatives  where 


?k  = £+  n 2 Ik’  k"1 Zk=l  ~n^~k  “ 2’ 


(5.12) 


it  has  the  coreesponding  non-central  distribution  with  the  non-centrality  parameter 


-1 


= q=1  Trace ( rR(  T(l:)  8 CR)  ), 


(5.13) 


where 


rk  -((  Yjk)  yjkJ.))  and  Ck  - limn^  nk^  , k=l,...,s,  (5.14) 

k 

which  we  assume  to  exist.  Asymptotic  optimality  results  hold  under  the  same 
setup  as  in  the  later  part  of  Section  4. 
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